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There is currently a great need for solid state lasers that emit in the infrared, as this is the op- 
erating wavelength regime for applications in telecommunications. Existing 7r-conjugated polymers 
all emit in the visible or ultraviolet, and whether or not 7r-conjugated polymers that emit in the 
infrared can be designed is an interesting challenge. On the one hand, the excited state ordering in 
trans-polyacetylene, the vr-conjugated polymer with relatively small optical gap, is not conducive 
to light emission because of electron-electron interaction effects. On the other hand, excited state 
ordering opposite to that in trans-polyacetylene is usually obtained by chemical modification that 
increases the effective bond-alternation, which in turn increases the optical gap. We develop a the- 
ory of electron correlation effects in a model vr-conjugated polymer that is obtained by replacing 
the hydrogen atoms of trans-polyacetylene with transverse conjugated groups, and show that the 
effective on-site correlation in this system is smaller than the bare correlation in the unsubstituted 
system. An optical gap in the infrared as well as excited state ordering conducive to light emission 
is thereby predicted upon similar structural modifications. 

PACS numbers: 42.70. Jk,71.20.Rv,71. 35.-y, 78.30. Jw 



I. INTRODUCTION 

TT-conjugated polymers have attracted wide attention 
in the past decade as light emitting organic semiconduc- 
tors. Light emitting diodes with 7r-conjugated polymers 
as emissive raaterials have reached very high efficiencies 
and lifetimeall, and there has also been considerable ac- 
tivity in the area of lasing using organic materials since 
the original demonstration of spectral na*ij«wing in these 
systems under intense optical pumpin^u'Q. One serious 
limitation in this area has been, however, that all light 
emitting 7r-conjugated polymers to date emit in the visi- 
ble or UV. Telecommunications use infrared radiation, so 
lasing at these wavelengths is desirableO. Whether or not 
TT-conjugated polymers can be synthesized that have light 
emission in the IR, and what should be the structural fea- 
tures of such a polymer are interesting theoretical ques- 
tions. Theoretical work on these topic&iwas recently initi- 
ated by one of us and his coUaboratorsBEl, where however, 
the emphasis was on explaining the unexpectedly high ex- 
perimentally observed photoluminesjcpice (PL) efficiency 
in polydiphenylacetylenes (PDPA)QoE3, which consist of 
a backbone carbon chain that is the same as in none- 
missive polyenes and side groups that are phenyl groups 
instead of atomic hydrogen. It was claimed that in an 
ideal long chain of PDPA the optical gap is smaller than 
that in trans-polyacetylene (t-PA). Real PDPAs emit iii 
the visible because of their finite conjugation laigthsO, 
but it was claimed in the earlier theoretical workoQ that 
relatively small optical gaps and strong PL would be ex- 
pected in any 7r-conjugated polymer where in addition to 
longitudinal conjugation there occurs conjugation over a 



few molecular bonds in the transverse direction. 

PDPA's, however, because of their complicated elec- 
tronic structures cannot be considered as the prototype of 
these novel class of nmterials. Furthermore, although the 
previous calculationsoLl were carried out to a high level of 
accuracy using the multiple reference doubles configura- 
tion interaction (MRDCI) technique, by themselves they 
do not give a clear mechanistic explanation of why the 
ordering of the excited energy states in PDPA's are op- 
posite to that in the linear polyenes such that light emis- 
sion becomes feasible In the present paper we therefore 
choose the simplest model system that shares the struc- 
tural feature of transverse conjugation with PDPA's, and 
analyse its lowest excitations. We demonstrate that finite 
transverse conjugation over a few bonds is the key to ob- 
taining light emission in the IR (provided, of course, long 
chain systems can be synthesized) . Our theoretical work 
shows that there exists a different mechanism for the en- 
ergy crossover between the lowest one- and two-photon 
states in 7r-conjugated polymers than the more-common 
one that depends on effective bond alternatiorilJlla. The 
simplicity of our model system allows us to perform nu- 
merical work iising a diagrammatic exciton basis valence 
bond methodli3 within which the basis functions retain 
their local character, and diagrammatic interpretations 
of wavefunctions can be given. A complete insight to 
the theory of excited state ordering in these complicated 
systems is thus obtained, which at the sajiae time fully 
supports the conclusions of our earlier workBEl that used a 
more traditional approach. To the best of our knowledge, 
although the exciton basis has often been used in the 
past for calculations on coupled two-level systems (often 
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though with severe approximations)" 
there e:> ii^t fev such calculations for coupled multilevel 



systemsc 



again sometimes within simpler models. 



Additional interest of our work then comes from its ap- 
plicability to realistic multilevel systems. 

Emission in the IR would require 7r-conjugated poly- 
mers whose optical gaps are smaller than or compara- 
ble to that of t-PA. Linear polyenes and t-PA are, how- 
ever, weakly emissive, because the lowest two-photon 
state,— the 2Ag, occurs below the optical 1B„ state in 
theseEj. The same_js also true for long-chain polydi- 
acetylenes (PDAs)B3. This is a conseqence of moder- 
ate electron-electron (e-e) interactions. h|Pt||vij'pa.-the tt- 
electrons, and is by now well understoodBElE3Ea-E£l. The 
optically pumped 1B„ in these systems decays to the 
2Ag in ultrafast times, and radiative transition from the 
2Ag to the ground state lAg is forbidden. Strong PL 
in systems like poly-paraphenylenevinylene (PPV) and 
poly-paraphenylene (PPP) implies excited state order- 
ing E(2Ag) > E(1B„) [where E(...) is the energy of the 
state], which is a consequence of enhanced bond alter- 
nation within the effective linear chain model for these 
systemstllE3. Within the effective model, the lowest ex- 
citations of systems containing phenyl rings as part of the 
conjugated main chain can be mapped onto those of a lin- 
ear chain with bond alternation that is much larger than 
that in t-PA. Explicit calculations then indicate that even 
with the same e-e interactions, chains with larger bond 
alternation can have E(2Ag) > E(1B„). Since enhanced 
bond alternation necessarily increases E(1B„), it appears 
that strong PL should then be limited to systems with 
optical gaps larger than that of t-PA, and emission in the 
IR from TT-conjugated polymers would be impossible. 

Our goal is to demonstrate that materials obtained 
by "site-substitution" of t-PA, in which the hydrogen 
atoms of t-PA are replaced with transverse conjugated 
groups will simultaneously have small optical gaps and 
E(2Ag) > E(1B„). The physical reasoning for this is 
as follows. Systems obtained by such site-substitution 
consist of units that are molecular (for example, trans- 
stilbene in the case of PDPA) , and consecutive units are 
linked by a single longitudinal bond. With moderate e-e 
interactions, the ground state can be thought of as cova- 
lent (all atomic sites singly occupied). Optical excitation 
involves intra- and inter-unit one-electron hoppings that 
generate double occupancies on the units. Since the unit 
consists of a large molecule in the case of the substituted 
polyene, this double occupancy occupies the antibond- 
ing molecular orbital of the unit, which is delocalized 
over the entire unit. Thus the effective on-site Coulomb 
interactions (the effective Hubbard interaction, hereafter 
Ueff) is smaller, and this can simultaneously give smaller 
E(1B„) (relative to t-PA) and E(2Ag) > E(1B„). This 
particular idea is related to a very similar idea in the 
area of organic conducting charge-transfer solids, where 
it has long been believed that by going to a larger organic 
molecule [for example from TMTTF, tetrathiafulvalene 
to BEDT-TTF, bis(ethylenedithio)fulvalene] the molec- 




FIG. 1: The model substituted polymer, poly- 
diethylenepolyacetylene. The numbers are indices for 
the atomic sites wthin each unit (see text). 



ular Hubbard U decreasesEj. The only difference in the 
present case is that we are comparing the bare Hubbard 
U (in t-PA) and a molecular Uef / (in the substituted ma- 
terial), and that the interunit one-electron hopping here 
is much larger (comparable to the intraunit hoppings) 
than in the charge-transfer solids, so that molecular ex- 
citon ideas have to applied with caution. 

We consider here the simplest possible 7r-conjugated 
polymer with conjugated sidegroup, the hypothetical sys- 
tem poly-diethylenepolyacetylene, shown in Fig. |l|, in 
which all hydrogen atoms of t-PA have been replaced by 
ethylene groups. The reason for choosing this hypothet- 
ical system is that precisely because of its simple elec- 
tronic structure, it can be considered as the prototype 
model for this class of (future) materials. In addition, 
in order to illustrate the theoretical result we are after 
we need a physical system for which, (a) accurate many- 
body computations can be done within a minimal basis 
(since for any realistic system with large molecular units 
full configuration interaction, FCI, with the complete ba- 
sis is out of the question) for relatively long chains, and 
(b) proper justifications can be given for ignoring basis 
functions beyond the minimal basis. Note that (b) is 
at least as important as (a), since it is known that the 
2A_pan occur above the 1B„ also because of incomplete 
CiSEil. Thus for our study to be complete it is essential 
that we prove that our result is not a consequence of our 
choosing the minimal basis. The symmetry and simplic- 
ity of our model system, along with the particular basis 
we have chosen allow us to achieve this. 

In the next section we give the Hamiltonian, parame- 
ters and the scope of our calculations. In section IH 



we 

discuss the complete exciton basis for oligomers of the 
system shown in Fig. |l], following which we discuss our 
minimal basis, which in this case consists of the highest 
occupied molecular orbital (HOMO) and lowest unoccu- 
pied molecular orbital (LUMO) of each unit. In section 
IV we present our numerical calculations for the unsubsti- 
tuted and substituted polyenes to show that there exists 
a broad range of Hubbard U over which the 2Ag is below 
the IBu in the unsubstituted polyene, but above it in 
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the substituted material of Fig. Wavefunction analfe 
sis within the exciton basis valence bond (VB) methodE^I 
also shows that the effective correlations in the substi- 
tuted polyene are smaller. Following this, in section 
we give very detailed heuristic but rigorous arguments 
which show that this result is not a consequence of in- 
complete CI. The mechanism of the 2Ag - 1B„ crossover 
in polyenes is by now understood; specifically, the corre- 
lated 2Ag state is dominated by a particular class of two 
electron - two hole (2e-2h) excited configurations, viz., 
covalent configurations that do not contain any double 
occupancy. Since the 1B„ consists of ionic configurations 
only, there must exist a U where the IB^ energy is higher. 
Because of the pictorial natures of the wave£iinctions in 
the diagrammatic exciton basis VB methodtj, it is pos- 
sible to classify precisely all possible one electron - one 
hole (le-lh) and 2e-2h exciton diagrams. As we show 
in section M, based on the symmetry properties of our 
simple model system this information can then be used 
to demonstrate that the higher energy MOs of the units 
play insignificant roles in the 2Ag - 1B„ crossover and 
that the minimal basis results are therefore accurate. 



occurs below the IBj, We perform similar cal- 

culations for oligomers of the the substituted materials, 
polydiethylenepolyacetylene (Fig. |^) here. For our hy- 
potheses in section | to be correct the following should be 
true. First, not only should the optical gap in the substi- 
tuted polyene be smaller than that of the corresponding 
polyene with the same number of backbone chain car- 
bon atoms at all J7, but the rate of increase of E(1B„) 
with U must also be smaller. The latter ia a definite 
signature of smaller effective correlation Ueff- Second, if 
Uef f is indeed smaller than the bare U in the substituted 
polyene, the bare Uc at which the 2Ag - IB^ crossover 
occurs in the substituted polyene must be larger than in 
the unsubstituted polyene. As a consequence, for a given 
e-e interaction it becomes possible to have the 2Ag above 
the 1B„ in the substituted polyene even as it is below the 
1B„ in the unsubstituted polyene. 



III. DIAGRAMMATIC EXCITON BASIS FOR 
POLYDIETHYLENE- ACETYLENE 



II. THE THEORETICAL MODEL 

We will consider our model polymer system within the 
simple dimerized Hubbard Hamiltonian, 

H^Hi, + ffee (la) 

Hle^- ty4^Cj- ^ (lb) 

-ffee = C^^f^i,T"-i,i (Ic) 
i 

In the above. Hie describes one-electron nearest neighbor 
hoppings of electrons, and Hee consists of the electron- 
electron (e-e) interaction within the Hubbard approxima- 
tion. The Fermion operator c| ^ creates a 7r-electron of 
spin a on site i, rii^^ is the number of electrons with spin 
tJ on site i, and (..) implies nearest neighbors. The one- 
electron hopping integrals are taken to be ti = 2.4(1 - 
d) eV, and t2 = 2.4(1 + 6) eV, with 6 = 0.07, correspond- 
ing to single and double bonds, respectively. These values 
of the hopping integrals are considered standard for tt- 
conjugated systems. We have ignored all Coulomb inter- 
actions other than the on-site repulsion U, since the 2Ag - 
1B„ crossover is related to this interaction only, with the 
spin-independent long range intersite Coulomb interac- 
tions merely modifying the magnitude of U at which the 
crossover occurs. Nonzero intersite Coulomb interactions 
make our exciton basis calculations much more compli- 
cated, without providing fundamentally new insight. 

Within the Hamiltonian of Eq. (1), the 2Ag occurs 
above the 1B„ at ?7 = in linear polyenes. For nonzero 
(5, as U is increased from zero, there exists a finite Uc at 
which the 2Ag - 1B„ crossover occurs, and the 2Ag now 



A. The complete basis 

It is clear that accurate many-body calculations for 
oligomers with the structure of Fig. |l| as the repeat unit 
can only be done within a reduced basis. The reduced ba- 
sis cannot be within the configuration space, as it would 
be simply incorrect to ignore atomic orbitals. A common 
approach in such cases is often to perform CI calcula- 
tions using a limited basis of molecular orbitals (MOs), 
but this is not suitable in the present case: MOs are 
delocalized over the entire space, and within the MO 
approach there is simply no information on the connec- 
tivity between the atoms that can allow choosing the 
reduced basis. It is for this reason that we have cho- 
sen the exciton basis VB method, which was previously 
used to perform j exac t calculations within the Pariser- 
Parr-Poplc modclE3H for linaar polyenes and to obtain 
diagrammatic wavefunctionsEj. As noted there, the ex- 
citon basis VB method is a hybrid of configuration space 
VB and the momentum-space MO pictures, and thus the 
basis functions within this approach retain the full in- 
formation on the connectivity between the atoms (as in 
configuration space) while at the same time, as in MO 
theory, the basis functions can be given certain hierarchy 
from energetic considerations. As we show later, both of 
these are essential for choosing the minimal basis as well 
as its justification. 

As in the case of polyeneJll we arrive at the exciton 
basis by rewriting the one-electron term Hi^ in Eq. 1(b) 
as a sum of two terms, one intra-unit and the other inter- 
unit, 

^ jjintra jointer 

The intra-unit part of the Hamiltonian is the simple 
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Hiickel Hamiltonian for the hexatriene units, 

5 

jj.ntra = -^,^1 ^^'1 + ^ + + h.c] (3) 

A' j = l 

where fi is an index for the units and j is an atom within 
the unit. The inter-unit part of Hi^ can be similarly 
written as, 

HIT' = -h E[^L4,.c,+i.3,. + h.c] (4) 

where we recognize that the bonds between the units 
involve the weaker hopping integral ti only, and only the 
atom 4 of the unit of the left unit and the atom 3 of the 
right unit are bonded (see Fig. |l| for the numbering of 
the atoms within each unit). 

The solutions to 77™*'"'^ are the intra-unit local Hiickel 
MOs, described by, 

6 

where A: = 1 - 6, with A: = 3 and 4 corresponding to 
the HOMO and LUMO, respectively. The reverse trans- 



formation, 

6 

fc=i 

where [Ajk] = [A^j]"^ gives the site-operators in terms 
of the local MO operators, and can be used to obtain 
Ijmter tcrms of thcsc, 

/J. k.k' ,(T 

(7) 

fjinter Operating on any exciton basis VB diagram 

leads to interunit charge-transfer (CT), which may in- 
volve pairs of bonding MOs, antibonding MOs, or bond- 
ing and antibonding MOs. The magnitude, and more 
importantly the sign, of the CT matrix element depend 
on the Ajk- In the case of linear polyenes, CT involving 
the bonding MO of the unit on the left lead to positive 
matrix elements of H™*^'^, whereas CT involving the an- 
tibonding-jMO of the left unit lead to negative matrix 
elementalj. The same is true here if k, k' in Eq. (|^) are 
limited to the HOMO and LUMO (fc, k' = 3, 4). 

Eq. (||) also enables us to rewrite the Hubbard inter- 
action as. 



M j=l kiM ks.ki 
M j ki,k2 k,ki,k2 

+ E I^J.feP •^i,fei^J\'=2^M,fe,iV'i,fci,T'/'M,fe2,T + E E ^J.'=i^J.fe2^J.'=3^i.fe4'/'|^,fei,TV'M,fe2,TV'|,,fe3,xV'M,A;4,i] (8) 

k,ki^k2 k\^k2 k^^k/^ 

I 



Here Nf^^kM = w^^k.a'^^i^.kM, the number of electrons 
with spin a in the fcth MO of the /ith unit. Thus Hp^ 
contains diagonal density-density terms, as well as off- 
diagonal density-electron transfer and electron transfer- 
electron transfer terms within the exciton basis. All 
terms, however, are local (i.e., intraunit), and electron- 
transfers due to the Hubbard interaction are strictly be- 
tween the MOs of the same unit. If the units possess 
both center of inversion and charge- conjugation symme- 
try (as is true here), an additional simplification occurs 
when calculations are carried out within a reduced basis 
containing-pairs of MOs related by charge-conjugation 
symmetryllj, viz., terms containing jproducts of density 
and electron transfer vanish exactlyli3. This is because 
the product Aj,kiAj^k2 i^i the density - electron transfer 
term has different signs for sites j that are related by 
spatial symmetry. Of course when considering the effects 



of the MOs not included in the reduced basis the effects 
of these terms have to be considered. 



B. The minimal basis 

The complete Hamiltonian, Hie + -ffee is way beyond 
our reach for any system containing more that two units 
of the substituted polyene of Fig. m At the same time, 
we show in section |V| that the two-unit oligomer is a 
special case, and is not of interest. We therefore work 
within a minimal basis which is chosen as follows. We 
first solve FJ^*''" and note the following: (i) the HOMO- 
LUMO gap is much smaller than the other bonding to an- 
tibonding energy gaps (see Fig. |2| for the single-particle 
energies), and (ii) Hl^*^^^ involves only atoms 3 and 4 
of each unit, and the contributions by the HOMO - 1 
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FIG. 2: The single-particle energies (in eV) for each isolated 
unit, with standard hopping integrals (see text). The dashed 
box indicates our minimal basis. 



and LUMO + 1 levels to the electron densities on these 
atoms {Aj^k in Eq. (||) for k = 2 and 5) are relatively 
small [\Aj[i\ = \Aj,e\ = 0.52, \Aj,2\ = \Aj,5\ = 0.21, 
\Aj^3\ ~ \Aj^4\ ~ 0.43, where j = 3 and 4], although 
the contributions by the outermost single-particle lev- 
els A: = 1 and 6 are large again (note, however, that 
the single-particle energy gaps between these last pair 
of MOs is rather large, from Fig. We therefore con- 
struct our minimal basis using only the HOMO and the 
LUMO and retain only k,k' ^ 3 and 4 in Eqs. ^ and |. 
Since these MOs are related by charge-conjugation sym- 
metry the density - electron transfer terms (but not the 
electron transfer - electron transfer term) in iJge vanish 
exactly within our minimal basis and the minimal exciton 
basis Hamiltoniands very similar to the one considered 
before for polyenealj, except that the Ajk are very differ- 
ent now (l^jfel = 2^^/^ in polyenes). Note in particular 
that with Hee hmited to the Hubbard interaction the re- 
maining electron transfer - electron transfer term in H^e 
consists of two electron transfers involving opposite spins 
only, and therefore virtual excitations involving the outer 
MOs do not occur. The lower occupied MOs therefore 
contribute a constant to the diagonal matrix elements of 
the exciton basis VB diagrams, which is the frozen core 
approximation. 

While single-particle energies and electron densities 
give guidance to the choice of the minimal basis, justifica- 
tion of the minimal basis requires the demonstration that 
the minimal basis can satisfactorily reproduce the results 
expected for the complete Hilbert space. This is easily 
done for the one-electron part of the Hubbard Hamil- 
tonian. Hie. We solve the full Hiickel Hamiltonian for 
oligomers with number of backbone carbon atoms N = 
4-10, and compare the exact lowest excitation energies 
with the approximate excitation energies obtained using 
the minimal exciton basis. Excellent agreement is found, 
as is shown in Table 1, where we compare the energies 
of le-lh B„ excitations for N = 4 - 10 (the le-lh Ag 
excitation energies can be obtained from the B„ energies 
and are therefore not given separately). 

Justification of the minimal basis even when Hee is 



nonzero is much more complicated. We postpone this un- 
til after we present our numerical demonstration within 
the minimal basis that the 1B„ optical gap is smaller in 
the substituted polyene even when U is included, and the 
Uc at which the 2Ag becomes lower than the IBu is much 
larger for the substituted system. As mentioned in sec- 
tion H, we are able to give a very complete proof that our 
reasult of higher Uc is not a consequence of incomplete 
CI, based on symmetry properties of exciton basis VB 
diagrams. 

C. The mechanism of the 2Ag — IBu crossover 

Although in reference |l^ detailed discussions of exci- 
ton basis VB diagrams for the case of linear polyenes 
were presented, in view of what follows in this and the 
next sections, it is useful to reemphasize some of these. 
We therefore discuss the mechanism of the 2Ag - 1B„ 
crossover within the the exciton basis VB approach. 
While it is known that the occurrence of the 2Ag be- 
low the IBu is a correlation effect that introduces strong 
configuration mixing between the "ground state" configu- 
ration and 2e-2h configurations, a key theme here is that 
only a specific class of 2e-2h states is actually relevant. 

Polyene exciton basis VB diagrams, and the substi- 
tuted polyene exciton basis VB diagrams within the min- 
imal basis are coupled by iJJ"*'^'" and Hee- We give here 
the complete expressions for Cu^s.o- and 0^,4.0- such that 
the efi:ects of Hl"*'^^ (see Eq. (m)) can be understood. 

C^,3,a = -0.52V'^,i^<, + 0.21V';.,2,<T + 0.43l/'^,3,a 

-)-0.43V'p,4,a - 0.21^^,5,^ - 0.521/^^,6,^ (9a) 
c^,4,tT = -0.52^^,1,^ - 0.21^^,2,(7 + 0. 437/^,3, 

-0.43^^,4,a - 0.21^^,5,a + 0.527/;^,6,. (9b) 

All numerical coefficients in the above are obtained 
from the invert of the Hiickel wavefunction matrix for 
linear hexatriene. The relative signs (though not magni- 
tudes) between the coefficients of the bonding and an- 
tibonding MOs in the site operators for the polyenes 
are the same as those between the coefficients of the 
HOMO and LUMO (^^,3,^ and ipf,,4.a) for the sub- 
stituted polyene in the above. Thus in the following 
both systems can be described simultaneously. Note 
that H"l*^"' from Eqs. ^ and (9) consists of terms 
k cr'^fj.+i-k' which Can be both positive or negative. 
These signs will be important in what follows in this sec- 
tion and in section 0. 

In Fig. 1^ we show the VB exciton diagrams and CI 
channels that are most relevant for the 2Ag - 1B„ 
crossover in a two-unit system, which can be either the 
two-unit polyene or the two-unit substituted polyene 
with the lower levels completely filled. Each arrow in 
this and similar figures below corresponds to one appli- 
cation of Hle*'^^ (Eq. (^) on the VB diagrams to the 
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States 


N = 


4 


N = 


6 


N = 


8 


N = 


10 


Exact Hiickel 


Exciton VB 


Exact Hiickel 


Exciton VB 


Exact Hiickel 


Exciton VB 


Exact Hiickel 


Exciton VB 


IBu 


1.81 


1.86 


1.41 


1.54 


1.17 


1.36 


1.02 


1.25 


2Bu 


3.37 


3.53 


2.48 


2.69 


2.22 


2.35 


1.99 


2.06 


3-Bu 






2.89 


2.83 


2.49 


2.37 


2.16 


2.08 


4B„ 






3.55 


3.72 


3.05 


3.18 


2.33 


2.66 


5B„ 










3.27 


3.22 


2.81 


2.84 


6Bu 










3.62 


3.34 


2.97 


2.91 



TABLE I: Exact versus minimal exciton basis energies (in eV) of the lowest le-lh Bu states in oligomers of the substituted 
polyene at U = 0. 




FIG. 3: CI channels and exciton basis VB diagrams relevant 
for the 2Ag - IBu crossover in a two-unit system with two 
MOs per unit. Bonding and antibonding MOs of each unit are 
occupied by 0, 1 and 2 electrons, with pairs of singly occupied 
sites coupled by singlet bonds (see text). Each arrow in (a) 
and (b) represents one application of H^^'^^ on the diagrams 
on the left, to obtain the diagrams on the right of the arrow. 
The plus and minus signs of the matrix elements are obtained 
from Eq. (9), and their magnitudes are not shown. The linear 
relationship (c) shows that the diagram with crossed bonds 
is a superposition of two other diagrams with the same or- 
bital occupancies but with noncrossing bonds. The diagram 
on the right in (d) is triplet-triplet (TT), with two triplets 
localized on the two units, (e) shows the overlap between the 
nonorthogonal 2e-2h VB diagrams. 



immediate left, to get the VB diagramis to the immedi- 
ate right of the arrow. As mentioned above, the plus and 
minus signs here are obtained fron the equations for the 
site-operators above. Aside from the signs the CI chan- 
nels should also include the magnitudes A^kA^k' , which 
can again be obtained from Eq. (||), but which are not 
very important for our discussions. 

The first diagram in CI channel (a) is the product 
wavefunction of the ground states of two noninteracting 
units, with all electrons in the bonding MOs, wtich we 
refer to hereafter as the Simpson ground statalJ. This 
diagram is symmetric with respect to the mirror plane 
between the two units, occurs only in the Ag subspace, 
and dominates the lAg wavefunction for small to moder- 
ate U . One application of 77™*'='" generates the "minus" 



linear combination of two charge-transfer (CT) diagrams 
(two because CT can occur in either direction). The rel- 
ative minus sign is obtained from Eq. (9). Here and in 
the following a singlet bond, drawn as a line connecting 
MOs k and k\ on the same or different units, is defined 

as '^'"\i^UA',u'A -_^UA',k'Am as heioA A 
second application of Hl"*^"^ , such that CT is now in the 
opposite direction, generates the 2e-2h charge-neutral di- 
agram with crossed bonds (the last diagram in (a)) which 
can again occur in the Ag space only. The "minus" sign 
in front of this diagram in the CI channel shown occurs in 
the case of both unsubstituted and substituted polyenes, 
and will be relevant for later discussions. 

Fig. |^(b) shows CI channels originating from the super- 
position of two Frenkel exciton diagrams with intra-unit 
excitations. The two Frenkel exciton diagrams are ab- 
sent in the Ag subspace, and the "plus" linear combina- 
tion shown occurs in the B„ susbspace. Operating with 
jjmter ^]^g Frenkel exciton diagrams generates the 
same CT diagrams that are generated in the second step 
in (a), but now with a relative plus sign between them. 
Because of this relative plus sign a second application 
of if™*'^'' (again promoting CT in the opposite direction 
from the first step) leads to an exact cancellation, such 
that the diagram with crossed bonds does not occur in 
the B„ subspace. 

The mechanism of the 2Ag - IB^ crossover is then 
as follows. The CT diagrams in Fig. |^ have unequal 
numbers of electrons on the units, and are hence nec- 
essarily ionic in the language of configuration space VB 
theory. The Frenkel exciton diagrams are also ionic, as 
can be ceafirmed by direct expansion of exciton basis 
operatorsll3. The diagram with the crossed bonds in 
Fig. §(a), however, is charge-neutral as far as the units 
are concerned, and as shown in Fig. ^c) is a super- 
position of two VB diagrams with noncrossing bonds, 
which can also be confirmed by expansion of exciton ba- 
sis operators. Equally importantly, a 2:1 superposition 
of these last two diagrams (see Fig. ||(d)) corresponds to 
two triplets, or triplet-triplet (TT) excitation on the two 
unitst^, which is covalent in the language of configuration 
space VB theory. As has been emphasized by previous 
authors, the 2Ag - IB,; crossover can then be understood 



1.5 - 

(a) 



) 1 


2 


3 4 5 
U(eV) 


6 


7 8 










IB/ 










2A 




















(b) 



'o 1 2 3 4 5 6 7 8 
U(eV) 



FIG. 4: The excitation energies E(1B„) and E(2A9) for the 
unsubstituted (solid Unes) and the substituted polyene (dot- 
ted hues), (a) N = 8, (b) N = 10. 



as a competition between CT and TT diagrams. Even 
though the diagram with crossed bonds is obtained in 
second order in Hl"*^"^ from the Simpson ground state, 
this second order process can be more important than the 
first order process leading to the CT diagrams, because 
of relationship (d) , since covalent diagrams have smaller 
diagonal matrix element of Hee than ionic diagrams. The 
2Ag is dominated by ionic diagrams at small U and oc- 
curs above the 1B„, while it is dominated by covalent 
TT diagrams at moderate to large U and occurs below 
the 1B„. We have discussed the crossover for a two- unit 
system only. The only new feature in long chains is the 
occurrence of CT and TT diagrams with long bonds be- 
tween units that are not neighbors (see reference |l3| and 
below). 



IV. NUMERICAL RESULTS WITHIN THE 
MINIMAL BASIS 

We performed exact numerical calculations for N = 4, 
6, 8 polyenes and quadruple-CI (QCI) calculations for 
the N = 10 polyene to determine Uc- Within the exci- 
ton basis the Hamiltonian matrix is nonsparse and this is 
the reason for doing QCI instead oiLexapt calculations for 
N = 10. As has been shown beforeaO, QCI is excellent 
for determining the energies of the lowest excitations, the 
IBu and the 2Ag, to high accuracy, at these N. Identi- 
cal calculations were done within the minimal basis for 
oligomers of polydiethylene-acetylene. In Fig. ^(a) and 
(b) we have plotted the excitation energies E(1B„) and 
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FIG. 5: The difference in Uc between substituted and unsub- 
stituted polyene as a function of 1/N. 



E(2Ag), with respect to the ground state energy (E(lAg) 
= 0) for N = 8 and 10, respectively. As claimed in sec- 
tion H, E(1B„) is lower for the substituted polyene in 
both N = 8 and 10, for all values of U. Note in par- 
ticular that E(1B„) for the substituted polyene increases 
with the bare U much less rapidly than for the unsub- 
stituted polyene, which is a distinct signature that Ueff 
in the substituted material is smaller. The Uc at which 
E(2Ag) becomes smaller than E(1B„) in the unsubsti- 
tuted polyenes is very slightly larger in N 10 than 
in N = 8, in agreement with previous worlS. Impor- 
tantly, the 2Ag remains above the IB-^ in the substituted 
polyene for a much larger range of the Hubbard U, and 
Uc for the substituted polyene is larger than that of the 
unsubstituted polyene by about 1 eV. The initial increase 
in E(2Ag) with U for substituted polyenes in Fig. |^ is real 
and has been seen previously in density matrix renormal- 
ization groupj-QDMRG) studies of long polyene chains by 
Shuai et. alEB (see, in particular. Fig. 1 in this refer- 
ence). As discussed there, this initial increase is related 
to the magnitude of the 2Ag - 1B„ gap at J7 = 0. The 
smaller this energy gap, the more the tendency to ini- 
tial increase in E(2Ag) (although E(2Ag) does begin to 
decrease again at large U). The substituted polyenes 
considered here have relatively small 2Ag - 1B„ gap al- 
ready at N = 10, which explains the behavior of E(2Ag) 
in Fig. ^. Taken together, the initial increase of E(2Ag), 
and the less rapid rise in E(1B„) give the larger Uc in 
the substituted polyenes. In Fig. || we have plotted the 
difference AUc between the Uc values for the substituted 
and the unsubstituted polyenes against the chain length 
N. AUc increases modestly between N = 4 and N = 8, 
but then is seen to have a tendency to saturation. This 
behavior of AUc is to be expected from what we have said 
in the above. The initial increase in AUc is due largely to 
the initial increase in E(2Ag) with U in the substituted 
polyene, which occurs for a larger range of U the longer 
the system is. This in turn is simply a consequence of 
smaller one-electron 2Ag - 1B„ gap as N increases. As 
the chain gets longer, the rate of decrease of the one- 
electron 2Ag - IBu gap itself decreases and the range of 
U over which E(2Ag) has initial increase in the substi- 
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tuted polyene tends to saturate, leading to a tendency of 
saturation in the Uc as well as AUc- Within the minimal 
basis therefore, AUc is at least as large as that at N = 
10 in the long chain limit. The larger Uc in the substi- 
tuted polyenes is a direct consequence of the delocalized 
HOMO and LUMO in the monomer, an effect which is 
included in the Ajk- 

Information on the effective Coulomb correlation can 
also be found from wavefunction analysis. In our previous 
work we had studied-linear polyenes within the Pariser- 
Parr-Pople modcEjS, for both small bond-alternation {5 
= 0.07, where the 2Ag is below the 1B„, and large bond- 
alternation {5 = 0.3), for which the 2Ag is above the 1B„. 
We had found that while the difference in the IBu wave- 
functions in the two cases was very small, the difference 
between the 2Ag wavefunctions was substantial: the 2Ag 
for small 6 is dominated by TT exciton basis diagrams, 
while for large S there occurs a strong admixing of le-lh 
CT diagrams. It is then expected that similar differences 
in the 2Ag wavefunctions would occur here for the un- 
substituted and substituted polyenes. In Fig. ^ we show 
the 1B,( wavefunctions for the N = 10 polyene and the 
corresponding substituted polyene for U = 6 eV, which 
is larger than the Uc for the polyene but smaller than 
the Uc for the substituted polyene. These wavefunctions 
appear to be very similar in nature. In contrast, the 2Ag 
wavefunctions for the two systems for the same U, are re- 
markably different, as seen in Fig. 0. While the polyene 
2Ag is predominantly TT and TT CT (diagrams which 
appear to be similar to the TT but in which the lengths 
of the two interunit bonds are unequal are referred to 
as TT (g) CTtS) there is much stronger admixing of CT 
diagrams in the case of the substituted polyene, indicat- 
ing clearly that the effective e-e correlation in the latter 
is smaller for the same bare U. CT diagrams are nec- 
essarilv ionic in the language of configuration space VB 
theoryll3£3, and the greater ionicity of the 2Ag wavefunc- 
tion in the substituted polyene is related to its initial rise 
in energy with UcB. 



V. BEYOND THE MINIMAL BASIS, THE ROLE 
OF HIGHER EXCITATIONS. 

A large body of previous work has established that in 
the case of linear polyenes low order CI using g-MP basis 
does not capture important correlation effectscj'Ell. Thus 
within the the PPP Hamiltonian, double-CI fails to find 
the 2Ag below the IBu for polyenes with N > 8, and 
either QCI or MRD-CI becomes necessary to obtain the 
correct excited state ordering. Similarly, fuU-CI calcula- 
tions within the Hubbard model for N = 12 using the 
six innermost bonding and antibonding-MOs also fails to 
find the correct excited state ordering.Ell Thus the very 
fact that the minimal exciton basis gives a 2Ag - 1B„ 
crossover at all is already significant. Nevertheless, since 
for each N in polydiethylene-acetylene we ignore more 
MOs than we retain, it might appear that the numerical 



result that Uc for this system is much larger than that 
of polyenes is a simple consequence of our approxima- 
tion. We show in this section that this is not so, and the 
higher exciton basis VB diagrams play an insignificant 
role in the 2Ag ~ 1B„ crossover, especially in long chains. 
There is a fundamental difference between the CI in- 
volving completely delocalized MOs and the CI involving 
the exciton basis VB diagrams. The latter possess local 
character and are more suitable for calculations involv- 
ing local Coulomb interactions. We emphasize that even 
though our discussions below are heuristic and based on 
physical arguments rather than numerical calculations, 
they are complete and rigorous. 

Our demonstration of the weak role of the higher ex- 
citon basis VB diagrams will be done in several steps. 
First, we look at simple polyenes more closely, to deter- 
mine that only a certain class of higher excitations can 
lower E(2Ag) with respect to the IBu- Based on this we 
arrive at two specific criteria for relevance of higher exci- 
tations in polydiethylene-acetylene, viz., that they should 
also involve crossed bonds (in other words, that the rele- 
vant configurations are those that are related to the high 
energy TT diagrams); and that they are coupled to the 
lower energy TT diagrams. This allows us to severely 
limit our search of relevant high energy configurations. 
We then go through the complete list of such configura- 
tions, both intra- and interunit, asymmetric and symmet- 
ric with respect to the mirror-plane symmetry. Because 
the two-unit oligomer is a special case (see below) we 
discuss this separately. Following this we discuss longer 
oligomers. Our basic conclusion is that the longer the 
oligomer, the more irrelevant are the high energy con- 
figurations. The discussion that follows is unavoidably 
complex, and we have therefore given a separate sum- 
mary at the end of this section. The reader uninterested 
in the details of CI theory and the exciton VB basis can 
directly go to this summary. 

Before we present our discussion of higher MOs in 
polydiethylene-acetylene, we present additional discus- 
sions on the 2Ag - 1B„ crossover in polyenes that go 
beyond Fig. ^. Double-CI is not suitable for correct eval- 
uation of E(2Ag) since quadruple excitations that are 
coupled to the doubly excited configurations, and that 
can therefore lower E(2Ag), are excluded in double-CI. 
The CI channel describing this is shown in Fig. ||(a), 
where the quadruply excited diagram for the two-unit 
case is obtained from two applications of H^^*^"^ on the 
VB diagram with crossed bonds in Fig. ||, once again such 
that the two CT processes are in the opposite direction. 
E(2Ao) is also lowered by going to larger N. As shown in 
Fig. 11(b) this is due to CI between the two diagrams with 
short nearest neighbor and long next-nearest neighbor 
crossed bonds, again through two applications of H™*'^^ 
in opposite directions. Very similar CI process occurs be- 
tween the final diagram in Fig. ^(b) and the diagram with 
even longer crossed bonds and so on. Two things will be 
relevant in our discussions below. First, even though di- 
agrams with long crossed bonds are not directly coupled 
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FIG. 6: The 1B„ wavefunctions for the N = 10 (a) polyene, and (b) substituted polyene within the minimal basis calculation, 
for [/ = 6 eV. Notice that the wavefunctions are similar in nature. 
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FIG. 7: The 2Ag wavefunctions for the N = 10 (a) polyene, and (b) substituted polyene within the minimal basis calculation, 
for U = 6 eV. Pairs of diagrams within the square parentheses occur as nearly 2:1 linear combinations, reflecting TT character. 
Note that the two wavefunctions are now very different, with CT contribution in (b) much larger than in (a). 




FIG. 8: CI channels that lower the enerey of the 2Ag. Ar- 
rows have the same meaning as in Fig. H. (a) The CI pro- 
cess between the two-unit symmetric 2e-2h excitation and 
the quadruple excitation, (b) Similar CI process in the three- 
unit case that gives the symmetric 2e-2h diagram with long 
bonds. 



to the Simpson ground state, their relative weights (more 
precisely, the relative weights of the TT diagrams with 
long bonds) in the exact 2Ag can be large (see Fig. |^). 
This is because of the degeneracy of all 2e-2h VB di- 
agrams with crossed bonds within Hee + Hl"*^°' which 
makes CI channels of the type in Fig. ^(b) extremely im- 
portant. Second, the relative signs of the diagram with 
nearest neighbor crossed bonds and the diagram with 
long crossed bonds in Fig. |^ are the same. The same sign 
between them signifies constructive interference between 
the CI channels in Fig. |^(a) and Fig. ||(b), and as seen 



in Fig. 1^, all TT diagrams have the same sign in the 2Ag 
wavefunction (TT Cg) CT diagrams, defined in the above 
as TT-type diagrams that are accompanied by CT, as is 
true for the pair within the first parentheses in Fig. 
can have opposite signs due to the additional CT), while 
the relative sign of the Simpson ground state is opposite 
(the relative sign of the Simpson ground state and the 
diagrams with crossed bonds are the same in the lAg, as 
suggested by the CI channels; see reference |l3| ). 

In order to understand the role of additional MOs in 
polydiethylene-acetylene we will use the underlying in- 
formation in Fig. ||. Two classes of 2e-2h excitations are 
relevant for the 2Ag - IBu crossover in the substituted 
polyene, intra- and interunit. The former are reached by 
application of Hee (all terms in Eq. (p^)) to the Simpson 
ground state, the latter once again through two appli- 
cations of H™^'^'^ . The intramolecular 2e-2h excitations 
that are relevant for the 2Ag ~ IBu crossover are known 
from the existing literature oh-CI calculations using the 
MO basis for small moleculeaHa. As far as the interunit 
high energy 2e-2h excitations are considered, the foUow- 
ii^ is obvious from the previous sections and reference 

The only interunit 2e-2h VB diagrams that are rel- 
evant for the 2Ag - IB^ crossover are the ones with 
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crossed bonds between units, since these are the only in- 
terunit diagrams that are charge-neutral. Strong CI in- 
volving diagrams with crossed bonds is essential for hav- 
ing contributions from higher energy TT diagrams. 

Having identified the class of 2e-2h diagrams witliin 
the complete basis that may be relevant for the 2Ag - 
IBu crossover we note the following. Unless the higher 
energy 2e-2h diagrams are directly coupled to the lower 
energy 2e-2h diagrams by ifge or iJJ"*'^'' they have a very 
weak effect on E(2Ag). In the absence of CI between the 
high and low energy diagrams, we have a nearly block- 
diagonal Hamiltonian, with each block containing its own 
set of excitations, with only the matrix elements of the 
Simpson ground state linking the blocks. Thus the 2Ag 
wavefunction continues to have strong CT contribution 
as shown in Fig. |. Thus E(2Ag) is lowered only if there 
exists CI channels analogous to those in Fig ^ that couple 
the high and low energy VB diagrams with crossed bonds. 
This leads to our second conclusion: 

Unless the high energy relevant 2e-2h diagrams, con- 
taining intra- or interunit two-excitations, are coupled to 
the low energy diagrams with crossed bonds through two 
applications of Hl"*'^^ , they play an insignificant role in 
the 2Ag - IBu crossover. 

We have limited ourselves only to second order pro- 
cesses between excitations (fourth order from the ground 
state) in the above, since if the higher energy diagrams 
are coupled in even higher order, their effects on E(2Ag) 
are even weaker than quadruple excitations involving the 
minimal basis. Also, the couplings have to be necessarily 
through iJJ"*'^'' , since all couplings through H^e are intra- 
unit only. The above restrictions allow us to severely 
restrict our search for relevant high energy 2e-2h excita- 
tions and still reach the correct overall conclusion con- 
cerning the role of MOs beyond the minimal basis. In 
the following we discuss the case of the two-unit system 
first, followed by the more general case of arbitrary chain 
length. Our overall conclusion is that because of the local 
nature of iJJ"*'^'' within the VB exciton basis, the effect 
of the couplings with the high energy 2e-2h diagrams in 
long chains is very weak. Because of the very detailed 
nature of the discussions that follow, we have given a 
summary of these discussions in the last subsection. 



A. Two units: Coupling to intra-unit 2e— 2h states 

Exact 2Ag wavefunctions for small polyenes, obtained 
by FCI usipff the MO basis have been obtained by Schul- 
ten et. al.Ej. The dominant components of the 2Ag in 
molecular hexatriene are shown in Fig. ^(a) - (d). Of 
these, (a) is a superposition of le-lh excitations that 
can only increase E(2Ac,). As shown in Fig. ||(e) this le- 
lh configuration is coupled to intermolecular CT states 
in first order by H™*'^^ . The 2e-2h configuration (b) 
has been retained in our minimal basis calculation. This 
leaves the molecular configurations (c) and (d) in Fig. ||, 
and as shown in Fig. ^ (f) and (g), these are coupled 




FIG. 9: (a) - (d): Intraunit excitations that dominate the 
molecular in hexatriene. All 6 MOs are shown now. The 
lowest three in each unit are bonding MOs, while the upper 
three are antibonding MOs. (e): CI channel between one of 
the single excitations of (a) and a minimal basis CT diagram, 
for the two-unit system, (f) and (g): CI channels between in- 
traunit 2e-2h excitations and the two components of the TT 
diagram within the minimal basis, also for the two-unit case, 
(h) CI channels between higher energy B„ Frenkel configura- 
tion and the minimal basis Frenkel configuration. Arrows 
have the usual meaning in (e) - (h). 



to components of the TT diagrams in the minimal basis 
by TJJ"*'^'' in second order. Thus both le-lh and 2e- 
2h higher energy molecular Kg excitations are coupled 
to the minimal basis Ag diagrams, and without actual 
calculation, it is difficult to say which of the two CI pro- 
cesses involving the molecular A^ states shown in Fig. ^ 
will dominate near U = Uc- Furthermore, there occurs 
also CI between higher molecular B„ configurations and 
the minimal basis Frenkel configurations, as is shown in 
Fig. |^(h), and this process will lower E(1B„), which will 
have a larger bandwidth upon including all MOs. This 
last result was also pointed out in references from 
transition dipole moment calculations. 



The most important conclusion that one draws from 
Fig. H is that all interactions between minimal basis dia- 
grams and the molecular diagrams are strictly local. The 
very nature of the CI process is such that the configura- 
tions that occur in the molecular 2Ag have matrix ele- 
ments only with minimal basis TT diagrams in which the 
two triplet excitations occur on nearest neighbor units. 
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(a): 



(c): 




FIG. 10: (a) and (b): Superpositions of asymmetric 2e-2h ex- 
citon basis VB diagrams with crossed bonds, for the two-unit 
case. There are many such diagrams and their superpositions 
occur in both Ag and B„ subspaces, as is shown for one of the 
superpositions in (c) and (d) (see text). 






FIG. 11: Symmetric 2e-2h VB diagrams with crossed bonds 
in the two-unit ohgomer of the substituted polyene, (a) is 
the minimal basis diagram, (b) - (f): high energy diagrams 
ignored in the minimal basis calculation. 



B. Two units: Coupling to intermolecular 2e— 2h 
states 



Two-unit VB diagrams with crossed bonds can be clas- 
sified as symmetric or asymmetric with respect to the 
mirror plane between the units. Both even and odd su- 
perpositions of the asymmetric VB diagrams are then 
possible, and these VB diagrams therefore occur in both 
Ag and B„ subspaces. Examples of such asymmetric VB 
diagrams involving MOs beyond the HOMO and LUMO 
are shown in Fig. |l^(a) and (b). The simplest way to 
confirm that the superpositions shown do occur in both 
subspaces is to apply iJJ"*'^'' twice on either the Simp- 
son ground state (for occurrence in the Ag subspace) 
or the Frenkel exciton state (for occurrence in the B^ 
subspace). The application of HJ^*-'^'^ in this multi-MO 
case is achieved through Eqs. and ^, and the signs 
originate from the different signs for each "01 fc crV''2,fc'.cr 
^2 k iy^^,k' -ly- This is done for the superposition of 
Fig. |ro|'(a) in Fig. [l^(c) and (d). The occurrence of the 
asymmetric VB diagrams in both Ag and B„ subspaces 
indicates that they play no significant role in the 2Ag - 
1B„ crossover, and including them in our numerical cal- 
culations would not have altered the value of ?7c in N = 
4. 

More important in the present context are therefore 
the symmetric TT diagrams. We have shown in Fig. ^ 
the corresponding diagrams with crossed bonds. These 
symmetric diagrams occur only in the Ag subspace and 
can in principle lower the energy of the 2Ag. VB diagram 
|ll|(a) is the lowest energy diagram with crossed bonds, 
and we are interested in determining whether this di- 
agram is coupled to any of the diagrams (b) - (f) in 
Fig. [l|. 

Three out of the five sets shown in Fig. |ll], sets (d), 
(e) and (f), cannot have matrix elements of iJ™*'^'' with 
VB diagram (a) in second order, from orbital occupan- 
cies alone. This is because two applications of H^*-'^'^ can 



change the occupancies of only two MOs, while the VB 
diagrams in (d), (e) and (f) differ from diagram (a) by 
occupancies of four orbitals. This leaves only the dia- 
grams in (b) and (c), which are coupled to the diagram 
in (a) once nonorthogonality of the VB basis functions 
is taken into consideration (see Fig. |(e)), as is shown in 
the two CI channels in Fig. |l|(a) and (b). The final VB 
diagram in Fig. ^ has nonzero overlap with diagram (a) 
in Fig. |ll|, and hence it might be thought that the CI 
processes shown can reduce E(2Ag). 

This is where the signs of the final diagrams in all CI 
channels become important. The minus sign in front of 
the final VB diagram in Fig. E3(a) and (b), taken to- 
gether with the overlap of -1/2 between this diagram 
with noncrossing bonds and diagram [Tl|( a), indicates that 
the diagram |ll|(a) is generated with a plus sign from these 
CI processes. This is in contrast to the CI process in 
Fig. ||(a) where diagram |ll|(a) is generated with a mi- 
nus sign. There is thus destructive interference between 
the CI channels and the overall absolute "yield" of the 
low energy diagram with crossed bonds, starting from 
the Simpson ground state, decreases when CI channels 
involving the two higher energy symmetric diagrams (c) 
and (e) in Fig. ^ are included. Conversely, the rela- 
tive weight of the Simpson ground state decreases in the 
2Ag, and E(2Ag) increases. This interpretation of the 
CI process is unconventional but is perfectly valid (as 
can be readily ascertained by constructing and solving a 
Hamiltonian matrix that includes all the diagrams and CI 
channels shown in Figs |^ and |ll|, to determine the extent 
of mixing between the Simpson ground state and the low 
energy diagram with crossed bonds when the interfering 
channels of Fig. |l^ are included) since the Hamiltonian 
matrix in the exciton VB basis is constructed csscntally 
from the linear relationships shown in Figs. ^ and ^ 

There is an alternate way to arrive at the same conclu- 
sion, which uses the relationship between VB diagrams in 
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FIG. 12: CI processes involving two of the four symmetric 
VB diagrams of Fig. ^ that are coupled to the minimal basis 
symmetric 2e-2h diagram of Fig. |ll|(a) (see text). The double 
arrows denote two consecutive applications of _ffj"*'^'", promot- 
ing interunit CT in opposite directions. Plus and minus signs 
are consequences of applying Eq. ^ Similar CI processes in- 
volving the other two symmetric VB diagrams of Fig. ^ also 
give the same final diagram with minus sign. 




FIG. 13: The absence of coupling in second order between 
high energy long bonded 2e-2h VB diagram with crossed 
bonds and minimal basis VB diagrams with crossed bonds. 
Two applications of Hil*"^^ on the diagram on the left does 
not generate either of the two diagrams on the right. 



C. Longer chains: Coupling to interunit 2e— 2h 

states 



exciton basis and configuration space. We have already 
remarked that the singly excited units in the Frenkel ex- 
citon diagrams of Fig. |3| are ionic. In the case of polyenes, 
each singly excited unit is the superposition of ionic con- 
figurations 20 - 02, where the numbers denote atomic 
site occupancies. Thus two singly excited units contain 
(20 - 02) (20 - 02), i.e., the ionic configurations 2020, 
2002, 0220 and 0202, and it is only the superposition in 
Fig. ^(d) that cancels these ionic contributions and gives 
the covalent VB diagram. Configuration mixing of the 
type in Fig. ^ therefore involves ionic diagrams at en- 
ergy 2U (as opposed to covalent diagrams at zero energy), 
and hence can only increase E(2Ag). 



Our overall conclusion then is that already at the level 
of two units, we find that interunit high energy dia- 
grams with crossed bonds play a weak role, and only 
couplings with intraunit excitations can perhaps lower 
E(2Ag) [note that such couplings will also lower E(1B„), 
see Fig. ^(h)]. In such a short system, the role of in- 
traunit 2e-2h excitations can be strong. If all MOs are 
included in a FCI calculation, it is possible that the true 
2Ag is below the true 1B„ in N = 4, although the dif- 
ference in their energies should be considerably smaller 
than in butadiene. Conversely though, if the molecular 
unit is such that the molecular 2Ag is higher in energy 
than the molecular IBu, our discussion indicates that the 
2Ag in the coupled system must necessarily be higher in 
energy than the corresponding 1B„. This is in agreement 
with our previous calculations on oligomers of PDPAQ'El 
that used a multiple reference double- CI approach using 
a completely delocalized MO basis, and that found that 
the oligomer 2Ag is above the oligomer 1B„. The unit 
in PDPA is trans-stilbene, in which the molecular 2Ag is 
higher than the molecular 1B„. 



For longer chains we focus on interunit 2e-2h VB dia- 
grams only, since all interactions involving molecular 2e- 
2h states are restricted to minimal basis diagrams with 
nearest neighbor TT bonds. As discussed in the above, 
CI processes of the type shown in Fig. ^(b) are par- 
ticularly important, because of the degeneracy between 
the short and long-bonded crossed diagrams. The nature 
of the exact 2Ag here and in reference |l^ both indicate 
this. One can then visualize such CI processes involving 
all pairs of MOs fc. A:', which generate diagrams with long 
crossed bonds in longer chains, with the relative weights 
of the diagrams with short crossed bonds spanning two 
units gradually decreasing with chain length. Whether 
or not these diagrams with long crossed bonds are cou- 
pled to the minimal basis TT diagrams then decide how 
strongly the calculated minimal basis Uc is affected. 

The classification of intermolecular 2e-2h VB diagrams 
as mirror plane symmetric or asymmetric persists in 
longer chains. Consider first the class of symmetric VB 
diagrams with long crossed bonds between distant units. 
Between any pair of molecular units there are as many 
of these as in Fig. 11, Since, however, if™*'^'' consists 
of nearest neighbor CT only, there are no symmetric di- 
agrams with long bonds that are coupled to the lowest 
energy diagram with crossed bonds in second order now. 
This is shown in Fig. |l^, where we have taken a three- unit 
version of one of the diagrams in Fig. ^ that did have 
coupling with the lowest energy diagram with crossed 
bonds. Thus the minimal requirement of having the same 
occupancy in two out of four MOs is no longer sufficient 
to give coupling to the low energy diagram. 

New CI channels become important now with the 
asymmetric diagrams with crossed bonds. A subset of 
the asymmetric diagrams can now give long-bonded min- 
imal basis TT diagrams, as shown in Figs [l4|(a) and (b), 
where we have given two such examples. It should be 
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FIG. 14: Examples of CI processes that can generate the 
long bonded minimal basis VB diagram with crossed bonds, 
starting from the Simpson ground state, with high energy 
asymmetric 2e-2h diagrams with crossed bonds as the inter- 
mediate states. The double arrows again denote two consec- 
utive applications of HJ"''^''. Note the opposite signs for the 
final product in the two channels, which is a consequence of 
charge-conjugation symmetry. 
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FIG. 15: Complete list of all asymmetric 2e-2h high energy 
VB diagrams with crossed bonds that are coupled to the 
Simpson ground state as well as the long bonded minimal ba- 
sis VB diagram of Fig. |l^, through CI processes of the type 
shown there. The plus or minus sign associated with each di- 
agram corresponds to the sign on the minimal basis diagram, 
when the CI channels as in Fig. ^ are constructed, starting 
from the Simpson ground state in each case. 



clear that the reverse of what is shown in this Figure, 
viz., generation of short-bonded minimal basis TT dia- 
grams from long-bonded asymmetric diagrams is equally 
probable. In general high energy asymmetric diagrams 
with bonds spanning L units can give low energy sym- 
metric diagrams with crossed bonds spanning L - 1 or L 
+ 1 units. We discuss only the former, since the essential 
physics of all these processes are identical, and originate 
from charge-conjugation symmetry alone. Since we are 
interested in two consecutive CT processes in opposite 
directions only, the particular subset of asymmetric di- 
agrams with short crossed bonds that can give the long 
bonded minimal basis TT diagram must satisfy a very 
strict condition, viz., both the HOMO and the LUMO of 
one of the molecules in the initial state must be singly 
occupied. For a reason that will be clear below, we have 
given in Fig. |l^ the complete CT steps, starting from 
the Simpson ground state, that show how the asymmet- 
ric VB diagrams are generated in second order, and then 
another second order process gives the minimal basis VB 
diagram with long crossed bonds. The CI between the 
high energy asymmetric 2e-2h diagrams and the minimal 
basis diagram with crossed bonds clearly seems to suggest 
that the former are relevant, according to the two criteria 
of relevance we have proposed, and can lower E(2Ag). 

There are two reasons why CI processes such as those 
in Figs. lU (and all similar processes in which asymmetric 
diagrams with crossed bonds of length L generate sym- 
metric diagrams with minimal basis crossed bonds with 
lengths L-l-1 or L-1) contribute very weakly to the ex- 
act 2Ag. First, the intermediate VB diagram with short 
crossed bonds in Fig. ||(b) is degenerate with the fi- 
nal TT state (as already pointed out in the above, this 
is what makes the contribution of long bonded TT dia- 
grams in the 2Ag of polyenes large), while both the inter- 



mediate 2e-2h states in Fig. 14 have diagonal energies 
that are different from that of the final TT state. Even 
more importantly, we note that while the "yield" of the 
long bonded crossed diagram in Fig. ^(a) is positive. 



that in Fig. |l|(b) is negative. The opposite signs re- 
sult from straightforward applications of Eq. and are 
a consequence entirely of charge-conjugation symmetry, 
which determines the signs in these equations. This has 
a very important implication for the extent of configura- 
tion mixing between the Simpson ground state and the 
long bonded TT diagram, viz., the configuration mixing 
induced by any one channel is reduced by the other. The 
simplest way to see this is to imagine that the two super- 
positions of asymmetric diagrams in Fig. |l^(a) and (b) 
have the same diagonal matrix elements of H^e + HYl*"^"" ■ 
In that case, by operating with H^f"^^ we would have sim- 
ply obtained the superposition of all four diagrams (the 
relative signs between which are determined entirely by 
Eq. (||)), and now a second application of iJJ"*'^'" would 
simply give zero. 

In Fig. |l^ we have given all the asymmetric diagrams 
with crossed bonds in the three-unit chain that can be 
obtained by applying Hl"*'^^ twice to the Simpson ground 
state, and that in turn give the long-bonded minimal ba- 
sis diagram with crossed bonds when further operated 
on by Hl^*'^^ (in each case we have given one diagram 
only, but superpositions of the type in Fig. 
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are be- 
ing implied). The positive and negative sign associated 
with each diagram is the sign on the long bonded final 
diagram when CI channels as in Fig. 14 are constructed 
with these diagrams as intermediate states. Not surpris- 
ingly, the total number of such diagrams in Fig. |l5| is 
even. Taken together with the fact that the configura- 
tion mixing between the Simpson ground state and the 
long-bonded TT diagram due to any one of these chan- 
nels is small to begin with (due to the nondegeneracy 
between the intermediate asymmetric state and the final 
state), the even number of channels, with half of them 
giving positive matrix elements and the other half giving 
negative implies a nearly total cancellation effect. Since 
the above result is a consequence of charge-conjugation 
symmetry alone, it remains valid for all sets of asym- 
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metric diagrams with crossed bonds, independent of the 
bond lengths. We therefore conclude that exactly as the 
symmetric diagrams with long crossed bonds, the asym- 
metric diagrams with crossed bonds have very weak CI 
with the minimal basis TT diagrams in long chains. 



D. Summary 

Our basic conclusion is that high energy intraunit 2e- 
2h excitations that have been ignored in our minimal ba- 
sis calculations can affect the final results. However, cou- 
plings between these and minimal basis TT diagrams are 
strictly local, and moreover, due to the nondegeneracy 
between the molecular 2e-2h states and the minimal ba- 
sis 2e-2h states the effect of this coupling on long chains 
should be weak. Recall, for example, that the relative 
weights of minimal basis TT diagrams with nonnearest 
neighbor bonds can continue to be relatively large even 
as the lengths of the intcrunit bonds keep increasing (see 
Fig. 1^ and also reference 13), but the couplings in Fig. ^ 
are only to TT diagrams with nearest neighbor bonds. In 
contrast to the intraunit 2e-2h states, the coupling of the 
minimal basis TT diagrams to the interunit high energy 
2e-2h VB diagrams is extremely weak. Indeed this lat- 
ter coupling can even raise E(2Ag) slightly in the case of 
two units. The weak coupling with high energy interunit 
states is a symmetry effect that will persist in any system 
which has spatial and charge-conjugation symmetries. 



VI. CONCLUSION 

In conclusion, we have analyzed the excited state or- 
dering in finite chain analogs of the simplest model conju- 
gated polymer with transverse 7r-conjugation within the 
dimerized Hubbard Hamiltonian. The system considered 
consists of a backbone polyene chain with ethylenic side- 
groups. We developed a multilevel exciton basis descrip- 
tion for the system, and within a minimal basis, which 
consists of the HOMO and the LUMO of each unit, we 
showed that the effective electron correlations in this sub- 
stituted polyene is smaller than that in the unsubstituted 
polyene. The reason for this is that each carbon atom 
of the polyene is now replaced with multiple atoms over 
which a double occupancy is delocalized, thereby giving it 
a transverse bandwidth that lowers its energyO. Smaller 
effective correlations, in turn, suggests that the 2Ag is 
above the 1B„ in this class of materials, and that this 
provides an interesting route to obtain organic polymers 
that emit in the IR. 

The chief advantage of the exciton basis approach is 
that the role of the high energy configurations that are 
ignored within the minimal basis calculations can be de- 
termined precisely, because of the local character of the 
basis functions, as well as the local nature of the con- 
figuration interaction within the Hubbard Hamiltonian. 
Our very detailed analysis here indicates that few of the 



high energy TT configurations that we ignore have CI 
in second order with the minimal basis TT configura- 
tions. In particular, the CI between interunit high and 
low energy TT diagrams is very weak. The only rele- 
vant CI between high and low energy states that can at 
all lower the energy of the polymer 2Ag involves the in- 
traunit configurations that describe the molecular 2Ag 
state. Note that this already ensures the following re- 
sult, viz., if the molecular unit is such that the molecular 
2Ag is higher than the molecular IB^j, then since within 
the minimal basis the 2Ag is also higher than the 1B„, 
the overall E(2Ag) is also higher than the overall 1B„. 
This in turn predicts that that E(2Ag) must necessarily 
be larger than E(lBt() in PDPA, since the molecular 2Ag 
is higher than the molecular IBu in the fundamental unit 
trans-stilbene. Our conclusion here is in agreement with 
the previous conclusion jajbout PDPA. which was arrived 
from computations onlyOtO. 



The molecular 2Ag configurations interact only with 
the minimal basis TT diagrams in which the two triplets 
occupy nearest neighbor units. In a long chain system, 
the relative weight of such TT diagrams with nearest 
neighbor bonds is small, since there exist in long chains 
TT diagrams with all possible bond lengths, and their 
contributions are comparable. At the same time, in 
the isolated real hexatriene molecule, which forms the 
unit in our model system, the lowest one- and two- 
photon states occur at a very high energy of ~ 5 eV, 
whereas the polymeric IB^ and 2Ag occur below 1.6 
eV, which is E(1B„) for t-PA. Thus any interaction be- 
tween molecular and polymer states decreases progres- 
sively with increasing chain length, and we expect that 
even for our model system the minimal basis results 
are qualitatively correct. Importantly, of course, our 
model calculations are for the illustration of the prin- 
ciple only, and real systems in which the molecular 2Ag 
is higher than the molecular 1B„ can always be devel- 
oped with suitable modification of the molecular unit. 
The real experimental challenge here is to develop sys- 
tems in which steric interactions are small enough that 
oligomers with long conjugation lengths and small opti- 
cal gaps can be synthesized. Finally, we mention that we 
have not discussed here known 7r-conjugated poly metes 
with low band gaps, such as p£>ly(isothianaphthene)E£l 
and poly(isonaphthothiophene)EZl. It is of interest that 
in these systems too there occur transverse conjugation 
over a few bonds in addition to the conjugation along the 
polymer axis, although the nature of the transverse con- 
jugation is slightly different from that in our model sys- 
tem consisting of a substituted polyene. The electronic 
structures and the small optical gaps in these-|Systems 
have been discussed within one-electron theoryEa and fu- 
ture theoretical work within many-electron approaches is 
clearly of interest. 
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